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ABS TRACT 


A group of stations in the North American Arctic region have 
been analyzed for statistical determination of temperatures at 
mandatory-pressure levels. For each station the temperature at a 
key level, peculiar to that station, has been forced in at the 
first step and retained at each subsequent step, in the development 
of stepwise regression equations which predict temperatures at the 
mandatory levels. In general, eight-step predictions, in terms of 
inter-level thicknesses, were found to give optimum specification 
of the desired temperatures. The best estimate of the regional 
atmosphere, which is conditionally dependent upon the existence of 
an extreme 1% probability of the forcing-level temperature, is 


obtained with a high degree of confidence. 
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I. INTRODUCTION 


Although many factors must be considered in determining the 
thermal environment of airborne ordnance, the most dominant are the 
flight conditions and the atmospheric temperature profile. While air- 
craft and missile performance capabilities are well defined, there is 
a great deal of variability in the atmospheric temperature profile 
which may be traversed by the vehicle. Therefore many model atmo- 
spheres have been developed for meteorological purposes, but none is 
completely satisfactory for use in aerothermodynamic heating 
calculations. 

Supplemental standard atmospheres have been developed by 
COESA(1966) to provide information on spatial and temporal variability 
across selected latitude bands and by seasons. However the supplement 
standards are intended to be means with respect to the latitude bands 
and seasons, and are not representative of the extremes which must be 
considered for aerodynamic design purposes. Richard and Snelling (1971) 
proposed the adoption of a revised MIL-STD-210 hot and cold model atmo- 
spheres which represent global extreme values at successive 2 km eleva- 
tions. These extremes are deviations from the station mean at the 
elevation involved, which are exceeded at the 1%, 54 and 104 probability 
levels, and were derived through analysis of extreme data from global 
meteorological temperature maps yee ae constant pressure levels. 
More recently a second draft, MIL-STD-210B, has been published (1972). 
Both of these extreme sets of profiles pose serious design problems when 


changes of elevation are involved. 





The most important limitation of the two recently proposed 
extreme-atmosphere sets is that they cannot be realistically utilized 
in conjunction with flight profiles and trajectories which involve 
changes in altitude. The temperature extremes at differing altitudes 
are in most cases representative of different geographical locations 
rather than of conditions occurring in real time over a real location. 
Polar and tropical standard and extreme atmospheres have been developed 
for such applications but are still subject to the deficiencies noted 


above. 


TABLE 1. Proposed locations of the 14 world-wide extreme temperatures 
at indicated levels (AFCRL, 1970), 


Pressure 
altitude 1% Cold extreme 14 Warm extreme 
(ft) location ‘ewes ealieyal 


Oymyakon, USSR Insalah, Algeria 
Oymyakon, USSR Insalah, Algeria 


Hall Beach, NWT * Babylon, Iraq 
Resolute, NWI “ New Delhi, India 
Thule, Greenland* New Delhi, India 
Thule, Greenland* Alert, NWT 
Karachi, Pakistan Alert, NWT * 
Singapore Thule, Greenland* 





The object of this study was to develop the most probable verti- 
Gal profile, conditional with the existence of a 14 extreme at a’ given 
level, for Gach station set forth in Table 1 (but here’ limited to those 
Nagoya aimeastc com och station or Lables i tas beensassociated 
with a listed 14 extreme for the indicated month according to a pre- 
liminary draft of the MIL-STD-210B proposed atmosphere, which is 
graphically illustrated in rig. 1. The asterisk-marked stations all 


occur in the Canadian Arctic and are either extremely cold at the 


10. 
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indicated level in winter or extremely warm in the summer stratosphere. 
This subset of the stations in Table 1 was isolated from the full set 
because the data were Metis ASAE ile in conveniently summarized 
form from checked rawinsonde data. The financial assistance and kind 
cooperation of the Commander, Naval Weather Service is gratefully 
acknowledged in acquiring the data on magnetic tape from the National 
Envinsonuentaie Data oervaec. sASivael le MNOGEM Cato lingaes Or each Station 
and month of interest, four years (1967-70) of atmospheric soundings 
were processed to determine the temperature means, standard deviations 
and inter-~level correlations required for application of the multiple 
regression analysis procedure to be used in the determination of the 


most probable temperature profile at the 1% probability level. 





Eis PBATA=“PROGCESSENG 


For each of the stations marked with an asterisk in Table 1, four 
years of checked rawinsonde data, for the period 1967 through 1970, 
was provided by the Naval Weather Service Environmental Detachment, 
Ashville. Those stations listed under the cold extreme of Table 1 
were subject only to January data samples, and those under the warm 
extreme, only to July samples. The individual sounding samples were 
arranged in the sequential format indicated in Table 2 and recorded 
on magnetic tape. Height and temperature data for the indicated 
pressure levels were arranged Ade six sets in each of the first three 
rows and five sets in the fourth row. Note that the data is arranged 
at 50 mb increments from 1000 to 200 mb inclusively and at somewhat 
smaller increments thereafter. 


TABLE 2. Arrangement of sounding temperature and geopotential 
height data by pressure levels k=1,...,23 in each sounding. 


Index k = i 2 3 4 5 6 
level (mb) 1000 950 900 850 800 750 


Index 


level 


Index 
level 





With the data in the format of Table 2, it is convenient to 


determine the mean temperature T@,) and istandard deviation om for the 


Us) 





indicated station eee specified by Table 1. The temperature TPS) 
is to be used in a multiple regression analysis to determine the most 
probable vertical temperature profile at all mandatory-sounding levels 
associated with that "forcing-level temperature" T@,). henetect of 
such station data is further analyzed for determination of the above 
statistical parameters for those soundings having 10% and 14 extreme 
temperatures at the forcing level. The T,¢.10) and T,¢.01) extreme 
temperatures can then be forced into the derived full-sample regression 


procedure to obtain the corresponding extreme profiles. 


TABLE 3. Nominal temperature extremes at indicated stations and 
pressure levels at the 10% and 14 values of a normal distribution. 


January Cases 


Full Sample 
Forcing Sample a1. 7. Cae) 7a 


Station evel Ty T ,¢.10) T,¢.01) Size Extreme 















Hall Beach r7 6G) -35.6 -42.3 2uS 26 








Resolute . T (500) —47.1 edhe) 237 25) 








Thule T (300) oar) I) os 213 ale 










Thule T (200) -63.5 -70.6 213 18 


July Cases 


Thule ot (100) = 1.5 mateo Re) 236 23 





TE IL SHO) 






T (200) 


Table 3 lists the forcing level pressures and sample sizes with 


the 10% and 14 limiting temperatures for each station under considera- 


tion. Estimates of the temperatures Ts), at the probability éxtremes 


14 





of 10% and 14 have neon computed assuming a Gaussian distribution at 
each J-level. These estimates, T,¢.10) and T,¢.01) Of Table 3, are 
derived using the computed full sample-means T(P,) and standard devia- 
tions Oy of the 1967-70 samples together with the aul ee arora 


relationships 


if 
| 
1 


+ 
T,¢.10) 1.28176, 


(1) 


+3 | 
ee 


T,¢.01) = 2.32670, 
where the plus sign is to be used for a warm extreme and the minus sign 
for a cold extreme. 

Analysis of the data revealed an insufficient number of soundings 
occurring at the 1% extreme of T(P,) for a test of the statistical 
significance of these results. It was decided therefore, to give 
primary consideration to the 10% extreme T@,) soundings and to deduce 
the properties of the 1% subset as a limiting case of the 10% sample. 
Thus the same multiple regression treatment used to predict the TCP) 


mandatory-level temperatures for the full sample has been performed 


upon the "10% subset" containing T(P,) extremes. 


15 





III. THE REGRESSION METHOD 


A. THE STEPWISE REGRESSION PROCEDURE 

Many problems in research require the extensive analysis of large 
amounts of data, and the process for handling it should be made as 
automatic and rapid as possible. To this end, the Biomedical Computer 
Programs (Dixon, 1966) were developed. One of these programs, BIMED 
O2R, provides stepwise regression analysis capability and will be 
utilized here. 

BIMED O2R computes, in a stepwise manner, a sequence of multiple 
linear regression equations. That single variable is added to the 
equation at each step which produces the greatest reduction in the pre- 
vious unexplained variance. This variable is also the one which has 
the highest partial correlation with the dependent variable at the 
particular step in the analysis of the variance. Equivalently it is 
the variable which would give the highest F-statistic upon entry at 
that step. The F-statistic upon entry ais is expressed at step Yr as 


(Dixon 1966) 


p Cjn-7) t% (C.8.V.,t) — % (C.5.V.,r-1) 
% (U.E.V. 50) 


where 


% (C.E.V.,r) is the percent cumulative explained 
variance at step r, 


% {C.&.V.,r-1) is the’ percent cumulative explained 
Variance at step r-l, 


% (U.E.V.,r) is the percent unexplained variance 
remaining at step r. 


16 





This study utilized a statistical model expressible in the form 
Ty = Cy af C,T(P,) ar CX (2) 
Here Ty in (2) is the temperature predictand at each of the eight 
mandatory levels chosen alternately as 


DSO eS 5 NS 


M DP To and T 


ve 2) 


but with Ty - TPs) s since the forcing-level temperature is used as the 


forced predictor for each of the eight mandatory level T,, specifications 


M 


for the station under consideration. Thus for each of the stations in 
Table 3, eight of the following nine dependent variables were specified 
for use in the multiple regression procedure: 


Ty» the temperature at 1000 mb, T(1000) 


T the temperature at 850 mb, T(850) 


9? 
Ty» the temperature at 700 mb, T(700) 
Ty» the temperature at 500 mb, T(500) 
Ty the temperature at 300 mb, T(300) 
‘Wes the temperatures at 200@mb, 1T(200) 


a the. t@mperature at 150 mb, 1(150) 


Ys 


T the temperature at 100 mb, T(100) 


9? 
Ts the Gonpe naeniee ak, 7/O-mb, 2 @io). 

Ts) is forced into the specification of each of the eight possible 

predictands in the manner symbolized by Eq. (2). Xe co lled samc rue. 

ture function, represents a linear combination of up to seven inde- 

pendent variables, comprised of layer thickness-values and the derived 

coefficients which meet the requirements of the stepwise least squares 


technique. The thickness variables Zi were generated within the BIMED 


O2R program by utilization of the so-called transgeneration technique 


17 





which takes the difference of successive tabulated heights in a 
sequential manner from the surface to 70 mb. This procedure produced 
twenty-two independent thickness variables which are directly related 
to the relevant mean virtual temperature in the layer (Pi Pray? by 
the hypsometric equation. These thickness predictors proved to be 
highly useful in the specification of Ty 

In the specification of the dependent variable Te a statement is 
permitted as to the maximum number of independent variables to be 
tested for admission to the regression equation. In the tests con- 
Sidered here, an upper limit of eight independent variables was found 
in general to give optimal specification of each of the dependent 
variables. Still another convenient feature of the BIMED O0O2R program 
is that the forced entry into the regression of the temperature value 
being tested as an extreme, the so-called forcing-level temperature, 
is allowed even though its contribution to the explained variance may 
become small at the final step in the regression analysis. 

Tables 4 through 10 list the variables selected and coefficients 
of Eq. (2) for each station extreme level considered. The listed 
structure functions i are often highly dependent upon the thickness 


variables nearly adjacent to the level represented by T However, 


M’ 


diversity of thickness-variable selection does occur. 
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TABLE 4, Full-sample regression equations, with associated thickness- 
variables at Hall Beach, NWT, for the specification of the January 
temperatures at indicated levels using T 5 (700) as the forcing-level 
temperature. 


Ty = ~268.4072 + 0.14737, ae la Xx) 

xy = 0.5952, 0.5322, = 0.5212, =P 0.0612, - 0.0352, . on 0.0812, 
-0.0592, . 

T, = -285.3345 - 0.0541T., + ale0 X, 

Xo = -0.014Z, - 0.0162, =F 0.3192, aR 0.3242, _ 0.0202. = 0.0082, 

T, = -293.3438 - 0.0755T, a cae ig 0) xX, 

xy, = 0.1062, ar 0.1092, 5 a5 0.101Z,, ar 0.0952, , - 0.0362, 4 

T, = -277.0920 + 0.0092T, + 1.0 7 

Xe = -0.0412,, + 0.1712, , =F 0.0942, . - 0.0142, ¢ 

Te = -244,0715 + 0.0905T., +91.0 Xe 

Xe = -0.0262, - 0.01/72, - 0.0152,, a5 0.0592, ¢ ar 0.2012, J _ 0.0262, , 

Ts = -262.7598 + 0.0300T 1 lO Xo 

Xo = -0.0202, _ 0.0372, ., i 0.1452, 2 “2 0.1142, 4 - 0.0152,,, 

Te. =a 249 O0G 0.0065T, ano Xe 

Xe = -0.0032,, = 0.0262, , or 0.0992.) =F 0.1052, - 0.0472, 

Ty = -254.3952 + 0.04477, +e 0 Xo 

Xo = -0.0222, ~ 0.0082, - 0.0052). - 0.0292... ta 0. 3082Z.,, 

xy is specified in terms of heights z as 

x) = 0.555 (2969-24 QQ) +f 0.532 (290972959) - 0.521 (26.9729) 
+0.061 (25, ,-Zeq) - 0.035(2569-2399) + 0.0812, 99-2550) 
~0.059(z, 59-21 95). 
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TABLE 5. Full-sample regression equations, with associated thickness- 
variables at Resolute, NWI, for the specification of the January 
temperatures at indicated levels using T,G500) as the forcing-level 
temperature. 


m = —313.9878 — 0.0034T, + 1.0 X 


1 4 1 
Xy = 0.4102, a 0.4262, = 0.3922., f 0.2332, ss 0.0102, , Se 0.0682, . 
-0.056Z,., 
T, = ~269.0266 + 0.02001, mee) Xy 
Ky = -0.0222, - 0.0092, aa 0.3252., 43 0.2842, = 0.0082, 
T. = -289.0706 = 0.0791T, gal ers X., 
XK. = -0.0902, 3 0.1692, + 0.1332, abe 0.1242, a 0.1702, - 0.0112, 
T. = =-262.3394 + 0.04217, = ioe Be. Xe 
Xe = -0.007Z, - 0.0112,, - 0.0342, 4 1% 0.1592, , az 0.1002, . - 0.0172, ¢ 
Te = -272.8584 + 0.02321), tole. O Xe 
Xe = -0.0112, . aL 0.0642, , az 0.2062, , = 0.0352, 
T, = -264.4431 + 0.01847, a el 8) Xo 
Xo = -0.0112,, = 0.0322, 5 aL 0.1332, at 0.1222, , - 0.0152, 
To. = -270.4817 + 0.01147, me or 8) Xe 
Xo = -90.0052, BY 0.0202, ay 0.0832,, We 0.1212, _ 0.0582, 
Ty = -269.8740 - 0.0072T, tee) Xg 
Xo = -0.0052,, az 0.0042, , = 0.0052, = 0.0402, a 0.3252,5 
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TABLE 6. Full-sample regression equations, with associated thickness- 
variables at Thule, Greenland, for the specification of the January 
temperatures at indicated levels using T (300) as the forcing-level 


temperature. 

lets? 2 Loe O0550T. 1b .0)K 

X, = 0.5342, + 0.5542, - 0.3642, - 0.037Z, - 0.0412, - 0.0222, , 
+0.043Z, . 

T, = -272.4780 + 0.00457, + 1.0 X, 

X, = -0.019Z, + 0.338Z, + 0.3152, - 0.0202, - 0.0302, 

T, = 272.8140 + 0.0093T, + 1.0 X, 

K, = 0.057Z, — 0.1182, + 0.1662, + 0.1622, + 0.108Z, + 0.1022, 

Wate 220 o2o) te OVOSEST eX, 

Ky GOV AZ a mOmO7aZr On MLO Ze oO d 252) = 0.04220, = On 0192, 7 

T, = -246.8040 + 0.0773T, + 1.0 X, 

Pee 00s ee O02 O27 et 0.0592). + 0.198Z5 2 = 00132. 

T, = -270.2495 + 0.00837, + 1.0 X, 

K, = -0.009Z, + 0.006Z, - 0.038Z,, + 0.156Z,, + 0.0812), 

ty = OIG) ES ay OCOtaESa Ss IO 

Ky = -0.006Z, + 0.0122, + 0.061Z,, + 0.114Z,, - 0.0362, 

Ty = -283.4604 - 0.0250T, + 1.0 X, 

Ky = -0.0062, + 0.0112, 5 OM Oz ZO O25 05277255 


Ze 





TABLE 7. Full-sample regression equations, with associated thickness- 
variables at Thule, Greenland, for the specification of the January 
temperatures at indicated levels using T , (200) as the forcing-level 


temperature, 

T, = -331.3064 - 0.26591, + 1.0 X, 

K, = 0.5262, + 0.5432, - 0.3462, - 0.0352, - 0.0392), + 0.0862), 
0.1562. 

T, = -272.2056 - 0.00381, + 1.0 X, 

X, = -0.019Z, + 0.337Z, + 0.3142, - 0.020Z, - 0.0292, 

T, = -273.0007 + 0.0049T, + 1.0 X, 

X, = 0.057Z, - 0.1182, + 0.1652, + 0.1622, + 0.1092, + 0.1022, 

T, = ~269.6765 + 0.0154T, + 1.0 X, 

X, = 0.0742, + 0.0742, + 0.118Z,, + 0.1242), - 0.0452, - 0.0062,, 

T, = -264.6812 + 0.05081, + 1.0 X, 

X, = 0.0092, - 0.055Z,, + 0.179Z,, + 0.091Z,, - 0.0192, , 

T, = -195.1660 + 0.2897T, + 1.0 X, 

X, = 0.0102), - 0.0212), - 0.091Z,, + 0.1612), + 0.0812), 

E, =§-283.8765 — 0.02627, 4+ 1.0 X, 

Xe = 0.0062, + 0.008Z,, + 0.060Z,, + 0.1142, - 0.0362.,, 

Ty = 280.3738 - 0.0026T, + 1.0 X, 

X, = 0.0092, + 0.0102, , = 0.02625, * 0.27725, 


No 
No 








TABLE 8. Full-sample regression equations, with associated thickness-— 
variables at Alert, NWI, for the specification of the July tempera- 
tures at indicated levels using T (200) as the forcing-level 
temperature. 


wy 


aT 


—273aGoll — 0, 32037, Fal0 xX 


6 1 

eres Ors O02 tO),0332,— 0.0272, + 0.028.) 
OOS, 

-304.2915 - 0.0173T, + 1.0 Xy 

0.0262, + 0.024Z, + 0.1422, + 0.3992, - 0.0092, - 0.007Z,, 
+0.0362Z, 
-273.5454 - 0.00347, + 1.0 X, 

aE eee Pepe ORC eee O 27 a Oe tz 4 Ol VLAZ 5 
-0.0072Z, 
-278.9707 + 0.03401, + 1.0 X, 

2 2 = 

COO eee ON OeceueCOS22) Onl 20z7e eG, 112z,,/— 0.05320 
0.0102, , 

-248,0954 + 0.1898T, + 1.0 xX, 

= Ome Pehaeeet ONO 2a 00322 5.1742). 4) 0.2642, 0) 0.0822, . 
-0.0402, 

-258.6167 + 0.04951, + Oe 

0.0152, = 0-010Z,, = 9.0042), ~ 0.047Z,_ + 0.155Z,, + 0.101Z,, 
- 0.0142, 

-272.3201 + 0.0112T, + 1.0 X, 

=0.006Z, + 0.0102, - 0.004Z,, - 0.004Z,, + 0.061Z,, + 0.134Z,, 
-0.063Z,, 

-273.5842 + 0.01381, + 1.0 X, 

Sor IO cue On ONOZ seai0 0222 0,0052,/ 4 0.0082... = 0.0132,, 
OO Alan 
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TABLE 9. Full-sample regression equation, with associated thickness- 
variables at Alert, NWT, for the specification of the July tempera- 
tures at indicated levels using T, (150) as the forcing-level 
temperature. 


T, = 207.5527 - 0.08987, + 1.0 X, 

KX, = 0.4242, + 0.3742, - 0.2482, + 0.0242, - 0.0452, ~ 0.0352, 
+0.0222,, 

T, = -329.3796 - 0.1059T, + 1.0 X, 

XK, = 0.0262, + 0.0242, + 0.1492, + 0.3872, - 0.008Z,, + 0.0552. 

T, = -274.1506 - 0.01047, + 1.0 X, 

Byp- 9. 0332, — 0.12724 +10. 1602, + 0.1682, + 0.1252, + 0.1192, 

T, = -271.5181 + 0.03107, + 1.0 X, 

X, = 0.0802 + 0.0842,, + 0.1122,, + 0.104Z,, - 0.0372,, 

T, = -276.9827 + 0.11627, + 1.0 X, 

KX, = 0,047Z,, + 0.032Z,, - 0.1752, + 0.2642,, + 0.081Z,, ~ 0.0252, 

T, = -259.7761 + 0.1023T, + 1.0 x, 

X_ = 0.0072, + 0.0382, + 0.1842,, -0.023z,, + 0.0162, 

Ty = ~285.4023 - 0.0316T, + 1.0 X, 

X_ = 0.0052, + 0.0102, - 0.0022,, + 0.0622, + 0.136Z,, - 0.0622, 

Ty = -273.5064 + 0.01357, + 1a X, 

Xy = -0.0072, + 0.0132, - 0.0042, + 0.007Z,, - 0.012Z,, + 0.2702,, 





TABLE 10. Full-sample regression equation, with associated thickness-— 
variables at Thule, Greenland, for the specification of the July 
temperatures at indicated levels using T, (00) as the forcing-—level 
ECnpC ra turer 


T = —215.4731 + 0.3238T. + 1.0 X 


i 8 Ih 

X) = 0.3602, SF 0.3742, = 0.0752. - 0.0342, _ 0.0412, a 0.0202, . 
-0.0212, . 

T, = i Oe 0.02727, a wees lee 6) Xo 

X, = 0.0472, + 0.0482, ot 0.1052, t 0.4012, = 0.0122, = 0.0042, 
+0.0182Z,, 

Ty = -319.8106 - 0.1370T, St tk X., 

X., = ~0.0812/, ats 0.1642, te 0.1562, + 0.1112, 25 0.1162, BE 0.0472,, 
-0.0312,, 

Ty = -276.6890 + 0.0313T, oe. OU Xx), 

X), = 0.0112, i 0.0922, a 0.0942, , aE 0.0962,, Si 0.0962, , _ 0.0332, 

T. = -266.3/04 + 0.07647, toa Xe 

X. = 0.0142, = 0.0732, . Sh 0.2142.) =F 0.0702, = 0.0162, . 

Te eer, 27 > 0.1/766T, mele) X 

Xe = 0.0192, , = 0.0352, ¢ + 0.0612), ai 0.1232, , ci 0.1032, . = 0.0442, 5 

Ts = -184.3870 + 0.0909T, Tet) a 

Xo = -0.0302, _ 0.0502, a. 0.0312), - 0.0142, . at 0.0682, _ tr 0.0082, , 
40.0212, 4 

Ty = Nee 0.020T, ve we) Xo 

Xo = 0.0112, =F 0.0022, - 0.0082, ., - 0.0032, , ae 0.2412,, 


25 








B. RELATED STATISTICAL PARAMETERS 
Available as outputs from the BIMED O2R program are several related 

statistical parameters such as: 

a. multiple Ro. after r steps 

b. istandard error of estimate, S.E. 

Cc. mean value, T 

d. ‘standard deviation, o 

e. F-statistic after r steps, F 

f. percent reduced variance, Ree 
The significance of Re im relation Fo S.E. si given by (Crow et allie 
1955) as 

(S.E.)* = 0° ((N-1)/(N-r-1)) (1-R7) 


where 


(1,-T)* / Ny. (3) 
1 


Q 
i 
ft Re 


i 
It should be noted that S.E. is a measure of the error in specification 
of TM after application of the regression equation estimator of 
moran 2). 


Lin ae is the variance of the temperature sample and 


N is the sample size 

i is the sample-element index 

r is the number of predictors selected 

Re. is the multiple correlation coefficient after r steps. It 
follows that the percent unexplained variance (1 - on for sample sizes 


of 213 to 244, as contained here, is closely approximated by 


(Say ay ae =i = ne 
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The percent explained variance rae also called the percent reduced 
variance, is then given by 


= Sa. Coan. 9 3)? (4) 


The mean value of Rae found in the prediction of the 1000 mb tempera- 
tures using the full-—data set, was 0.8478 and 0.7480 with the 10% 
extreme-data set. For all levels above 1000 mb, the pooled mean Ro for 
the full-data set was 0.96/72 and 0.9596 for the 10% extreme-data set. 
The scores at 1000 mb are somewhat lower than those of the upper levels 
and this is to be expected, in that the surface variabilities bring 
periodic and random-effect terms into the specification of the 


temperature at 1000 mb. 
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ives ANALYILTCAL PROCEDURE US ING) THE® REGRESSTON BECHNIOUE 


A. RESULTS OF THE REGRESSIONS 

Part (a) of Tables 11 through 17 lists the mean, standard deviation, 
multiple correlation SECC eter: and standard error of estimate for the 
full-data sample of each station and at all mandatory levels except the 
forcing level. These statistics were computed after application of the 
appropriate regression equations listed in Tables 4 through 10. Part 
(b) of Tables 11 through 17 gives similar results, after application of 
the identical regression equations used in part (a), to a regression 
analysis of the data set which contains only the "10% extreme" set of 
forcing—level sample atmospheres. This means that T(P,) was chosen if 
it lay within the particular 10% Gaussian distribution extreme under 
consideration. The remainder of the mandatory levels of the extreme 
sounding, T,(-10); were then subjected to the same regression analysis 
as the full-data samples of part (a). This nominal set of extreme 
atmospheres was selected by choosing those soundings for which 
TCP.) > T,(. 10) for the warm extremes, and T(P,) < T,(.10) folk stiewcold 
extremes. These extreme 10% values T,¢.10) are defined by the Gaussian 
ecient lope limitsml@sted imekq. (1). 

Part (c) lists the mean temperature profile defined at mandatory 
levels which correspond to the set of soundings having T(P,) at the 1% 
extreme. There were two ways of selecting the 1% extreme cases: (i) 
application of an a priori test to determine the values T(P,) which 
fell within the required range indicated by the Gaussian distribution; 
Sci aLethiere were meevalues of TP) Potting ties cOnstraame, then 


either the two or three coldest (for cold extremes) or warmest values, 
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were chosen a posteriori from the regression-determined 1% set as repre- 
senting the set of soundings corresponding to the 1% extreme T(P )- 
values. The corresponding set of TOP) values were then simply averaged 
and Per idered to represent the expected vertical profile having T(P 5) 


in the 1% extreme, 


B. THE TEST FOR SIGNIFICANT SOUNDING DIFFERENCES 

In view of the small sample sizes in parts (b) and (c) of 
Tables 11,..., 17, it was not feasible to test for significant differ- 
ence by random sampling of the 10% set relative to the 1% set of T(P 5) 
generated soundings. It was decided therefore to test for significant 
differences in defining the ie values predicted in each part (b) of 
Tables 11,..., 1/7 against those ie values of part (a) obtained by the 
identical prediction process, but using the full-data sample: that is, 
in case (a) there were no a priori limiting conditions set up regarding 
a particular value of TPs). 

In order to check the latter hypothesis, a Student's t-test was 
performed to test whether the 10% samples of prediction means of part (b) 


were significantly different from those of the full-data sample means 


bone a) eee poolcd statistic relating Costhic fe 


nM prediction 


differences of the part (a) and part (b) means may be formulated after 
Brooks and Carruthers (1953) with (N+tn-(rt+2)) degrees of freedom, after 
computing the two sets of predictive mean values and their standard 


Purors .welnemt-StatlLstie LS formulated as followse 


, L 
= (N- 9)o~ + (n- 2)0° 1 ] 
t(Nin-(r42)) = (F “10/ | E,.10\ (44 
FULL? T 10 sini cad ay (i n i (5) 





Here the T's are the part (a) and the part (b) sample temperature means 


De (at the mandatory levels) based upon identical prediction equations. 
o, is the part (a) standard error of estimate, and Op 19 is the part 
9 ¢ 


(b) standard error of estimate. The number of degrees of freedom asso- 


ciated with the two OF; values are (N-r-1) and (n-2), respectively, in 


parts (a) and (b) with 


N = the full-sample size 
n = the nominal 10% extreme-sample size 
r = the number of predictors, r = 8. 


The pooled t-statistic may be computed from the horizontally adjacent 
sections of Tables 11,..., 17 and with the proper number of (Ntm-(r+2)) 
degrees of freedom for t. Tables 18 and 19 of Appendix 1 list the 
resulting values at mandatory levels for each station under consideration. 


These t-values were tested against critical t ,,-values (e.¢g., 


O1 
Crow et al, 1960) with the same number of degrees of freedom used in 
(5). This number of degrees of freedom usually fell into the tabular 
class for t-values corresponding to an infinite number of degrees of 


freedom. Thus for example t 1? with an infinite number of degrees of 


0 
freedom, is t = 2.326. Tables 18 and 19 show one value below the 99% 
confidence limit specified by t o1 = 2+ Oconee aat valic. correspondims 
to the 1000 mb level of Thule (with Py = 300 mb), still yields a con- 
fidence limit of 92%. The random chance of the iy SHeciitedtl1ons of 
at cera andes of! tablest ley... 17 being drawn from the same popu- 
lation if t = 7.0 is approximately one chance in 5 X iors and is much 


Siiallmmrotmteny sOem—eriogt OL Lhest-valuece am TablaselS and 19 


exceeded t = 7.0 with exceptions principally at P = 1000 mb. 
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C. THE TYPE-(C) PROFILE CLASSIFICATION 

As noted previously, it was not feasible to employ a t-test utiliz- 
ing the extremely small sample size of part (c) of Tables 11,..., 17. 
eee it is useful to explain physically the pattern of trends in 
the sample means Ty, of part (c) relative to those of (b). The code 
symbol "D" for "direct", used in part (c), means chee lence is dis- 
placed directly in the same sense from T,,(b) as T,,(b) is from T,,(a) at 
the level under consideration. The code symbol '"M" for "mixed" means 
that parce) = F(a) | < IT, (b) = ie) so that a consistent displace- 
ment in the temperature extreme does not occur in progressing across the 
columns of T, values. This code symbolism was used for all cases con- 
Siaered and is) summarized, tinder part (c) of Tables 11,..., 17. To 
some extent, the nature of these displacements, whether direct or mixed, 


is indicated also by Fig. 2,..., Fig. 8, which show only the comparison 


Ot be ec) and T,, (a) profiles. 





ee” Sa. -- 


V. RESULTS SUMMARY FOR INDIVIDUAL STATIONS 


A. WINTER CASES 
This section will consider those cases investigated as cold ex- 
tremes and pennant Eee Cables SR. we Lomand Fios, 7a. 
1. Hall Beach, NWT; T(P,) — T, (700) 
Table 1l(a) shows the results of the stepwise regression 
applied to the specification of the January full-sample temperatures 


(Sig) Epp BRD Beene. 7 


M and Ty) at the mandatory levels of 


8 
Hall Beach by an equation of form (2). An eight-predictor equation for 


i at each level, was generated by use of the BIMED O2R program with 


M? 
the results shown. Very high multiple correlation coefficients were 
obtained, indicating a functional relationship for Ty at all levels. 
Table 11(b) lists analagous results for the data sample corresponding 
touthe 1a, cold extremerat the forcing level, with T,(b) being deter- 
Mined by the Same specification equation used in part (a). Part (c) 
presents the resultant estimated mandatory-level means T,,(c) derived 
from the five cases which fell into the 1% extreme forcing level sub- 
set of part (b). Part (c) then, gives the expected vertical profile 
corresponding to the 1% extreme at T(P;) = T(700). It is suggested 
that the mean T,,(c) profile, relative to that of = Oe should be phys- 
ically as well as statistically explainable from the conditions imposed 


by the choices of P At Hall Beach, the latter selection at /00 mb 


ny 
was consistent with the existence of an extreme-cold troposphere, 


LAG) < T(b); and an extreme-warm stratosphere for winter. The only 


os 


case of a mixed classification M, occurs for T., (500) ; where according 
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FIG, 2. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Hall Beach, NWI, corresponding 
to the set of 14 cold extreme occurrences of T,(700). The thin 


solid line depicts the January mean Hall Beach vertical sounding 
(1967-70). 
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to Fig. 2, the crossover between the mean and extreme soundings occurs. 
The nature of SMe) profile is then consistent with the displacement 
of the T, (b) = Te profiles except at the crossover point near the 
tropopause. 
2. Resolute, NWT; TP) = T (500) 

in the case of Resolute, Table 12(c) shows) a'somewhat different 
result from that discussed above for Hall Beach. At Resolute we now 
start with the premise of an initially specified very cold 
a) < T 91 009). This level is excluded from being an effective 
crossover point in the results of Fig. 3 except in the region of P=300 mb 
where a slight tendency toward a warm-stratosphere does occur. Compar- 
ing the two cases, the selection of an extreme cold T, (500) dictates the 
presence of a cold stratosphere coupled with a cold troposphere below. 
For nearly all levels at Resolute, chem 7p = ee trend Gears out 
that indicated sr a Ms) hy in the sense of maintaining the ae 
atmosphere relatively cold at virtually all levels. In contrast with 
the situation at Ilall Beach, the Resolute type (c) atmosphere more 
nearly resembles the case of a cold tropospheric vortex with an over- 
lying stratospheric cold vortex, whereas the Hall Beach type (c) case 


seems to occur with the warm stratospheric anticyclone aloft. 


oye thule, Creenland: T@,) = T, (300) 

This case is graphically depicted in Fig. 4, which shows that 
the a priori condition of a 14 cold temperature extreme at P =300 mb at 
Thule is associated with a cold stratosphere as well as a cold tropo- 
sphere at all levels above and below P = 300 mb. The extreme vertical 
profile for this regime TCP 5) = ~(300)) at Thule resembles that at 


Resolute, where there was also evidence of a similar coupling between 
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FIG. 3. The heavy solid line shows the mean regression-determined 


vertical temperature sounding over Resolute, NWT, corresponding to 


the set of 1% cold extreme occurrences of T CS00}e) The thin solid 


line depicts the January mean Resolute sounding (1967-70). 


42 









100 


150 


2 


2? 


200 


— Be 
= +} a 
a c 
ag “a 
~ |. Z 
o 
9 300 — 
i] 
6 
500 
5 
700 
850 
1000 
-100 -80 -60 -40 -20 0 20 40 
Temp (°C} 


FIG. 4. The heavy solid line shows the mean regression-determined 

vertical temperature sounding over Thule, Greenland, corresponding 

to the 1% cold extreme occurrences of T,(300). The thin solid line 
depicts the January mean Thule sounding (1967-70). 
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FIG. 5. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Thule, Greenland, corresponding 
to the 1% cold extreme occurrences of T,(200). The thin solid line 
depicts the January mean Thule sounding™ (1967-70). 
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cold stratospheric temperatures in case (c) and a cold-extreme upper 
troposphere. Now, however, Fig. 4 shows no crossover point between the 
mean and extreme profiles. Likewise Table 13 shows the direct transi- 
tion D-type classification applicable at all levels of the mean 
temperature profile T,,(c). 

4. Thule, Greenland; TP.) = T (200) 

This case differs from the previous case (T(P 5) = T,(300)), in 
that the latter was associated by regression development with T,(200) = 
-62.83°C, whereas the coldest 14 input temperature T, (200) Wace= 7070s 
The imposition of the latter condition into the regression analysis at 
Thule (T (200) ) is associated with the coldest stratospheric regime 
encountered during the entire study of this paper. With T, (200) = 
=F FO all regression-generated stratospheric temperatures are typi- 
cal of those in the cold Arctic vortex at and above 200 mb. However 
Fig. 5, shows that in this case there is a crossover in the T,(a), T,(c) 
profiles. Furthermore by comparison of Table 14(c) with Tables 1l(c), 
12(c) and 13(c), it is clear that this case corresponds to the warmest 
of the four winter tropospheres discussed here. 

5. Combined Winter Regimes, Summary 
The analysis of the four sets of type (c) winter extremes has 

indicated that the 300 mb level is a key indicator in the type of atmo- 
sphere to be expected in the Canadian Arctic region. In summary, 
these regimes have been identified by the four stratifications: 

(a) Hall Beach T (700) , cold extreme 

(b) Resolute T (500), cold extreme 

(c) Thule T , (300) , cold extreme 


(d) Thule T (200), cold extreme 
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Case (a) corresponds to an abnormally cold troposphere with a 
resultant abnormally warm stratosphere (Fig. 2). Cases (b) and (c) 
correspond jointly to cold stratospheres and cold tropospheres (Figs. 3 
and 4), but neither cold layer exists at a particular extreme. Finally 
case (d) is characterized by a record cold stratosphere and a near- 
record warm lower troposphere. Cases (a) and (d) show negative corre- 
lations of temperature across the tropopause, a property documented at 
mid-latitude stations by Cole and Nee (1965). However cases (b) and 
(c) have a consistent positive correlation throughout the atmosphere 


and have yet to be documented fully. 


Be cUMMER CASES 

This section will consider those cases as having warm extremes at 
forcing levels and are summarized in Tables 15, 16 and 17 and Figs. 6, 
7 and 8 The same stepwise regression procedure was applied to the 
July data of the warm extreme stations as was applied to the January 
data of the cold extreme stations. The results are presented below. 

ieee sert. NW; T(Py) = T (200) 

Table 15 lists the statistical results of the analysis of Alert 

with T(P,) = T, (200. Consistently high multiple correlations are 
shown indicating a high degree of confidence in the predicted 2, values 
including those of part (c). This confidence is reinforced by ae con- 
Sistent direct transition D-type classification applicable at most 
levels of the mean temperature profile Tc) Mixed tendencies occur 
at 100 mb and 70 mb but are characterized by extremely small tempera- 
ture differentials (less than 0.5°C) in the AG. ~ T(b) trends. A 


crossover of the T, Cc) extreme and T,, 6a) profiles exists between 500 mb 
M \ 
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The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Alert, NWI, corresponding to the 
i auarmmextreme occurrence of [_(200). The thin solid line depicts 
the July mean Alert sounding MCT 70e 
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and 300 mb (Fig. 6), indicating an extreme warm stratosphere overlying 
a relatively cool troposphere, This phenomenon is also characteristic 
of mid-latitude negative inter-level temperature correlations shown by 
Cole Pn Nee (1965) to exist near the level of the tropopause during 
pemlods OL cyctonc—anticycione variability. 

ye o.ert, Nw; T(P) = T (150) 

Table 16 corresponds to the analysis of Alert with T(P y) = 
T (150). Very high multiple correlation coefficients in both parts (a) 
and (b) are again shown yielding a resultant high degree of confidence 
in the diagnostic procedure of isolating "extreme" atmospheres. 

Figure 7 indicates that a marked similarity exists between the Ar (2) 
and Te) profile pairs and the corresponding pair for Alert of T, (200). 
lablewlosindicates that the Be) regression-generated atmosphere is 
displaced primarily in a "direct'' sense relative to the iA) atmo- 
sphere. It should also be noted that the forcing extreme T (150) at 
the nominal 10% and 1% values differ only by about 1°C. The comparison 
between jee) and Hh ©) is generally of the same character as that of 
Hee) Belative sto 7, (2), CE lea). 

The same conditions of warm stratosphere and relatively cool 
troposphere, as well as the crossover between 500 mb and 300 mb, sug- 
gest a similar acnimareskeor with regard to this climatic regime as with 
the preceding one (Alert T (200) ). 

3. Thule, Greenland; TP) = T, (100) 

The final table of regression statistics, Table 17, corresponds 
to the analysis of Thule at TCP») = T (100). Once again there are very 
high multiple correlation coefficients in both parts (a) and (b). In 


addition the atmospheric profiles, both in the mean and 12 extreme 
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FIG. 7. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Alert, NWT, corresponding to the 
1% warm extreme occurrence of T,(150). The thin solid line depicts 
the July mean Alert sounding (1567-70) . 
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FIG. 8. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Thule, Greenland, corresponding 
to the 1% warm extreme occurrence of T._.(100). The thin solid line 
depicts the July mean Thule sounding (1967-70). 








(Fig. 8), were nearly identical to the corresponding profile pairs in 
the previous two warm cases (Figs. 6 and 7). This indicates a strong 
tendency of the synoptic processes acting at Alert in July to also be 
acting at Thule during the same period. 
4. Combined Summer Regimes, Summary 

Each of the July cases is based upon the assumption of a warm 
extreme forcing-level temperature existing at or above 200 mb. In each 
case, there proved to be a concurrent warm temperature at 300 mb, that 
is, stratospheric warming above 200 mb has extended downward to 300 mb 
in July. At tropospheric levels, however, the mean temperatures of 
parts. (b) and (c) are generally colder than the full~sample means. 
Figures 6, 7 and 8 are graphical presentations of the mean temperatures 
shown in parts (a) and (c) of Tables 15, 16 and 17. 

Three factors should be recognized as being relevant to the 
July data considered. (i) In two of the three cases, the nominal. 1% 
extreme temperature at the forcing level T(P,) did not occur, and the 
class was redefined by a posteriori methods to comprise the set of two 
warmest observations at the forcing level. (ii) The difference between 
the parts (b) and (c) mean mandatory level profiles (T,,(b) and T,,(c)), 
at and above 300 mb, never exceeded ieOner (iii) The dominant feature 
of all of the summer case-studies was the consistent occurrence of 
stratospheric-tropospheric temperature reversals noted in the 


preceding paragraph of this section. 


D1] 





VI. CONCLUSIONS 


The multiple regression procedure presented for development of 
vertically consistent model atmospheres is statistically dependent upon 
the input condition of a known extreme temperature at a given flight 
level. The procedure has been shown to be highly definitive in deter- 
mining the most probable vertical temperature profile over the stations 
of interest. Extremely warm stratospheric temperatures are not asso- 
ciated with extremely warm tropospheric temperatures, thus integrated 
ballistic density anomalies tend to be nearly compensated. Likewise, 
extremely cold tropospheric temperatures tend to be reversed in the 
stratosphere. These two points become extremely important in the 
determination of ballistic peramecer sl As a result, aerothermodynamic 
calculations can be made for flight conditions involving changes in 
altitude by using an atmospheric temperature profile which is consis- 
tent with the proposed MIL-STD-210B value at the primary level of 


interest, but is more realistic and less severe at all other levels. 
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